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Abstract. We study several properties of blazars detected in the gamma-ray energy range by 
comparing the EGRET sources with a sample of radio blazars which can be considered possible 
gamma-ray candidates. We define three classes: non-gamma-ray blazars, blazars with quasi-steady 
gamma-ray emission, and gamma-ray blazars with substantial activity level. By combining the 
information of detected and candidate AGNs, we characterise the blazar activity, including the 
discovery of a region of consistency between the gamma-ray flaring duty-cycle and the recurrence 
time between flares. We also find a possible relation between the activity index of FSRQs and their 
black hole mass. 



INTRODUCTION 

Among Active Galactic Nuclei (AGN), blazars show strong flux variability at almost all 
frequencies of the spectral energy distribution (SED). The EGRET instrument on-board 
of CGRO detected, above 30 MeV, blazars as a class of y-ray sources O, identifying 
67 objects, and detecting 27 candidates. Gamma-ray blazars are characterised by high 
variability on different time-scales, from one day (e.g. PKS 1622-297) to one month 
(e.g. PKS 0208-512). This large spread in time variability and the sparse coverage 
obtained by EGRET make it difficult to quantify parameters such as the duty-cycle (i.e. 
the fraction of time spent in a flaring state) and/or the characteristic time-scale of y-ray 
activity for this class of sources. 

A preliminary result was obtained by |@], based on the Second EGRET Catalogue 
assuming that the link between the radio and the y-ray emission can be used to derive 
the y-ray luminosity function. Under these assumptions, they found that a duty-cycle 
£ = 0.03 is consistent with the data. One weakness of this analysis is that the radio- 
gamma relationship is highly uncertain. 

The aim of this Paper is to summarize a recent result on the estimate of the blazar 
duty-cycle and y-ray variability [9] . 



THE BLAZAR /-RAY VARIABILITY 



In order to investigate the blazar y-ray duty-cycle, we analysed the Third EGRET Cata- 
logue (3EG) AGN sample, looking for recurrent activity. For each source we considered 
the number of viewing periods (VPs) during which it was close enough to the pointing 
direction (< 40°, the FOV radius) that its exposure for that VP was greater than zero. 
For each source we calculated the Exposure (EXP), defined as the sum of the exposures 
during each single viewing period. 
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FIGURE 1. Four different examples of blazar y-ray activity as a function of exposure. Top-left panel: 
3EG J0845+7049 with HSN = 1 despite its high exposure; Top-right panel: 3EG J1733 — 1313 which, 
for a similar exposure, has HSN — 5; Bottom-left panel: 3EG J0530+1323 which showed a possible 
periodicity in its high-state status, and Bottom-right panel: 3EG J0340— 0201 which has a very low 
exposure (~ 64 x 10 7 cm 2 s), but showed 2 high-states in a few observations. Crosses show the value of 
the discrete \j/ for each VP containing the source. 

In order to estimate the y-ray activity level, we computed the number of times that 
a source was in a high-state (HSN) i.e. when the flux of the z'-th VP was greater than 
a given threshold. We establish that a source is in a high-state when 1) the flux of the 
z'-th VP is in excess of 1.5 times the mean flux of the source, and 2) the 1-c uncertainty 
of the measurement at the high state is less than the deviation of the measurement from 
the mean. The mean flux is defined as the weighted mean of all the detections and upper 
limits, assuming the weight w( o) = 1/(7, where o is the uncertainty of the measurement 
for detections, while o = UL for the upper limits. If two high-states were less than 2 
weeks apart (the typical duration of an EGRET pointing), we gathered them into a single 
high-state coincident with the first one. 

Given EXP and HSN, we computed an activity index (i/a) of each AGN as follows: 



V/ src = HSN src X [EXP SI 
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This index provides an estimate of blazar y-ray activity, weighted by the exposure. 
Sources with high y have a y-ray activity which is highly variable. Sources with i// = 
do not show significant flux variability. 

We analysed the behaviour of each single source having yr ^ during its time history. 
Fig. [l] shows the discrete i// for some EGRET sources as a function of EXP. We 
show four different examples of activity behaviour: 1) 3EG J0845+7049, with a single 
HSN = 1 despite its large exposure; 2) 3EG J1733 — 1313, which, for a similar exposure, 
has HSN = 5; 3) 3EG J0530+1323, showing a possible periodicity in its high-state 
status, and 4) 3EG J0340— 0201, which has a very low exposure (~ 64 x 10 7 cm 2 s), but 
showed 2 high-states in a few observations. 

Sources that have HSN = 1 (e.g., 3EG J0845+7049) show a rise of \^ rc followed 
by a decay after reaching a single high state, while sources that have high HSN (e.g. 
3EG J0530+1323) show a saw-tooth pattern, corresponding to alternately increas- 
ing/decreasing of \jQ c . 

We distinguish sources for which the activity level ty/ = (population A, 36 per cent), 
and the ones having i// ^ (population B, 64 per cent). A remarkable difference appears 
in the two populations when we compare the ratio between flat-spectrum radio quasars 
(FSRQs) and BL-Lac objects (BLs). Population A is almost equally divided between 
FSRQs and BLs (~60 and ~40 per cent, respectively), while Population B is dominated 
by FSRQs (~ 80 per cent). FSRQs are generally more variable and more luminous in 
the y-ray energy band than BLs, which could explain the higher fraction of these sources 
in Population B. A possible correlation is found between the activity level and the black 
hole mass, with a probability P ian d of a randomly distributed sample of objects on the 
order of ~ 10~ 2 . We suggest that the use of the \]/ statistic is more robust than the use of 
the HSN, since the latter is not weighted by the EGRET coverage. 



THE BLAZAR GAMMA-RAY DUTY-CYCLE 

We develop a simplified model of EGRET blazar activity that yields rough, but quan- 
tifiable, estimate of physical parameters of interest. In this simplified model, we make 
several assumptions: 

1. that all of the EGRET blazars exhibit the same basic behaviour. While observation- 
ally unjustified, this assumption is necessary to model the average behaviour. 

2. that only those blazars detected by EGRET are part of the population of interest. 

3. that the behaviour of all EGRET blazars can be characterised by a simple model 
with only two free parameters, the duty-cycle, %, and the characteristic time-scale, 
T. Each blazar will spend a period of time whose average length is T at a low flux 
level (off) before emitting a y-ray flare of duration x at a high flux level (on), then 
returning to a low level. The duty-cycle is then defined as the fraction of time spent 
in the on state: 

Note that all blazars have the same characteristic time-scale T, but the duration of 
each individual quiescent periods is drawn from a Poisson distribution with mean 



T. 



We explored the parameter space of (x, T). For each pair of % and T under consideration, 
we generated 100 sets of 67 simulated light curves covering the entire time interval 
of EGRET Cycle 1 to 4 (1620 days), one for each blazar. The time sequence of on- 
off states was determined by drawing the durations of the off states from a Poisson 
distribution with mean T. The durations of the on states were fixed at x. To obtain the 
observed fluxes, we used a bootstrap procedure. Each light curve was compared with the 
observation history of the corresponding AGN. If a given viewing period (VP) coincided 
with an on state, the observed flux was randomly drawn from the distribution of all fluxes 
detected from all EGRET blazars during a single VP Otherwise, the observed flux was 
randomly drawn from the distribution of all upper limits (i.e. non-detections) at EGRET 
blazar positions during a single VP. For each AGN the flux values assigned to each VP 
were used to calculate its HSN (and therefore its categorisation as Population A or B) in 
the same manner as the actual EGRET blazar HSN. The result is 100 sets of simulated 
blazar light curves, each with a ratio Rab = #Pop.A/ (#Pop.A + #Pop.B). Thus for each 
pair of values of (x, T), we have obtained a distribution of 100 simulated values of RaBi 
each with a mean and variance. Therefore, we are able to define a confidence region 
of (X,T) pairs whose mean values of Rab are within one standard deviation of the 
observed value Rab = 24/67 = 0.36. For values of Rab close to the observed value, 
the simulations give a standard deviation of 0.05. We ran 100 simulations each for 390 
{Xi T) pairs, ranging from < T < 450 days and from < % < 1. Fig. |2] shows the 
results. Points in Fig. |2] represent the mean value of Rab for each value of T and X- 
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FIGURE 2. Mean ratio R AB = #Pop.A/(#Pop.A + #Pop.B) from simulations of EGRET -detected 
blazars only. 

The filled triangles represent those ratio values between 0.31 < Rab < 0.41, small filled 
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circles represent 0.0 < Rab < 0.31, while crosses represent 0.41 < Rab < 1-0. Taking 
into account our simplifying assumptions, including the restriction that the observed 
EGRET blazars form the entire population of y-ray emitting blazars, the only region 
consistent with the observations is the region where 0.31 < Rab < 0.41. The region with 
high Rab produces too few flaring sources, while the region with low Rab produces too 
many. This implies that the true characteristic time-scale T lies between 50 and 400 
days, the duty-cycle % becoming more extreme for lower values of T. Note that for the 
Stecker & Salamon duty-cycle of 0.03 the characteristic time-scale T would be on the 
order of 50 days. 

If we should take into account a wider sample of y-ray blazar, by including radio 
blazars emitting at a level of 250 mJy, we would obtain that the low Rab region to remain 
excluded, because adding more sources is unlikely to reduce the absolute number of 
Population B sources predicted, which is higher than the number observed. On the other 
hand, the (%, T) pairs in the high Rab region may now be consistent with the EGRET 
data. Therefore, the region in (x, T) space where formerly the Rab was consistent with 
the observed value now forms a lower bound to the possible values of T for a given value 
of*. 



THE BLACK HOLE MASS - ACTIVITY RELATIONSHIP 

Recently, estimated the black hole mass for 36 EGRET AGNs from their optical 
continuum luminosity; they claim a positive correlation (significance level 96.6 per cent) 
between the maximum y-ray source luminosity and the black hole mass. 

We look for a possible correlation between the y-ray maximum luminosity, L™ ax , 
and the black hole mass, Mbh, obtained by adopting the black hole masses reported 
in ||8j] for all the EGRET blazars except for 3C 273, whose Mbh was taken from 
B5J. The correlation is absent when considering the whole sample (Spearman statistic: 
p = 0.44, and P = 0.15), and becomes marginal if we take into account the FSRQs 
only (p = 0.60, and P = 0.08). The relation between luminosity and black hole mass 
is, approximately, Ly Mgpj for FSRQs, consistent with the predicted relation by 01, 
although their assumptions (y-ray emission is isotropic radiation produced by the steady 
and Eddington-limited accretion without considering the beaming effect) do not agree 
with recent results on y-ray blazars [4]. When analysing the possible correlation between 
the activity level, i//, and the maximum y-ray luminosity, we find no correlation (p = 
0.14, and P = 0.39). We found a possible correlation between the activity index \y and 
Mbh- Fig. |3] shows our results based on ten blazars with non-zero ty/ and the black 
hole mass estimation. The Spearman correlation parameters turn out to be p = 0.76 
and P = 0.01 1. This possible correlation seems to indicate that more massive BHs can 
induce higher y-ray activity, although our sample is too small to draw a firm conclusion. 
More reliable BH mass estimates of a larger blazar sample will be crucial in verifying 
the existence of this correlation. 
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FIGURE 3. y-ray activity (i/a) vs. black hole mass (M Bli ) for EGRET Pop. B blazars. A marginal 
correlation is present (p = 0.76 and P = 0.01, Spearman statistic). The solid line represents the best fit 
line, log (!//■) = 0.013 x log(MBH) — 0.099. Horizontal bars represent the typical black hole mass estimate 
uncertainty as reported by Qj. 



CONCLUSIONS 

The results presented here summarize a more detailed analysis presented in [9] where 
we discuss the implications on the duty-cycle estimate of a sample of candidate y-ray 
AGNs extracted from recently published radio catalogues. A firm estimate of the duty- 
cycle will be a scientific goal of future y-ray missions such as AGILE and GLAST . Their 
wide FOV (~ 60°) will allow monitoring of a large number of AGNs for each pointing 
on long time-scales. For the time being, an extensive black hole mass estimation by 
means of optical spectroscopy on FSRQs and BL Lac objects will allow us to better 
understand the relation between the central engine and the y-ray activity. 
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